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solar energy and waste heat. In recent years, researchers proposed that ionic liquids (ILs) as novel 
alternative absorbent combined with refrigerant such as water, ammonia, alcohols, and hydrofluor- 
ocarbons can be used as working pairs for absorption refrigeration cycle, heat pump, and absorption 
Keywords: power cycle. In this paper, researches done in imidazolium IL working pairs regarding to status of 
Tonic liquid | evaluation and selection methods, thermophysical property measurement and modeling, as well as their 
Working pair future prospect assessments, i.e., developing potential studies about the absorption cycle performance 


Absorption cycle : : ; z 
Thermophysical property adopting new working pairs, have been reviewed. 
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1. Introduction 


1.1. Significance of R&D of absorption cycle and ionic liquid working 
pair 


To deal with the changes in global climate, rapid growth of 
energy demand and rapid development of renewable energy and 
energy utilization, much attention has been focused on the 
research and development (R&D) of absorption cycles, such as 
absorption heat pump and absorption refrigeration cycle, because 
they can make effective use of low-grade heat | 1-4]. However, the 
R&D of absorption technology faces a series of problems. For 
example, the temperature range of the heat source is wider, with 
a range of 80—90 °C collected by economic and efficient flat solar 
collectors to 500—600 °C exhausted from gas turbines or internal 
combustion engines. The functions and configurations of new 
absorption cycles are more diverse and complex, such as cogen- 
eration systems that combine power and refrigeration cycles [5-7], 
compression and absorption hybrid cycles driven by low-grade 
heat and electrical energy [8,9], and distributed energy systems 
that combine cooling, heating, and power cycles [10] involving 
energy storage [11-15] and absorption technologies [16-18]. The 
above developments present new demands for the R&D of work- 
ing pairs of absorption cycles as a key and essential work. 

Generally, the working pair of absorption cycle is a multi- 
species solution, which is usually a binary solution consisting of 
two species with a large boiling point difference. The species with 
relatively low boiling point, i.e., the volatile species, acts as 
refrigerant, whereas the species with high boiling point, i.e., the 
non-volatile species, acts as absorbent. The main types, existing 
problems, and requests for future development of the absorption 
cycle working pairs are briefly summarized and listed in Table 1. 
Based on the species of refrigerant, the working pairs are divided 
into the following five categories: water, ammonia, alcohol, 
halohydrocarbon and hydrocarbon (HC). The combinations of five 
kinds of refrigerants and different absorbents can basically meet 
the needs of different applications. For example, the H20/LiBr 
working pair is mainly applied in room air conditioning. The NH3/ 
H20 working pair can meet cooling requirements below 0 °C. The 
two working pairs described above are the most widely used 
systems with environmentally friendly features. 


Table 1 
R&D of absorption cycle working pairs. 


Refrigerant Absorbent Existing problems Main of the past R&D 
species species 
H20 LiBr Corrosion, crystallization Adding other auxiliary me 
octanol and other species 
NH3 H20 Toxicity, difficulty in separating NH3 from H20 Other organic ammoniates 
1) Absences of techno-economic 

g significance g : f 

HFC & HC Organic 2) Values of ODP? and GWP? cannot meet New refrigerants with higher ODP 
solvents the needs of developing and GWP values 
Alcohol (CH30H, Organic To be devoid of techno-economic significance Few 
C2H50H, TFE) solvents as like as the system H20/LiBr 


* Ozone depletion potential. 
> Global warming potential. 


However, the H20/LiBr system presents certain disadvantages, 
such as crystallization, corrosion, and negative pressure operation. 
Some researchers have attempted to improve the H20/LiBr work- 
ing pair by adding auxiliary species [19,20]. The NH3/H20 system 
also has some drawbacks, such as difficulty in separation and 
toxicity. Some researchers have added non-volatile species to 
reduce the separation difficulty [21-23]. 

Alcohol working pairs are mainly composed of low carbon 
alcohols, such as methanol and ethanol, or fluoroalcohols, such as 
2,2,2-trifluoroethanol (TFE) and hexafluoroisopropanol (HFIP), as 
refrigerant and salts or high-boiling point organics as absorbent. 
Alcohols are inflammable and methanol is toxic, so greater interest 
has been given to fluoroalcohols. Fluoroalcohols are non-corrosive, 
non-combustible, and have good thermal characteristics. TFE com- 
posed of working pairs with some high-boiling point organics have 
been used, such as tetraethylene glycol dimethyl ether (E181) [24,25], 
N-methyl-2-pyrrolidone (NMP) [26], N,N’-dimethylpropyleneura 
(DMPU) [27], N,N’-dimethylethyleneurea (DMEU) [27], and quinoline 
[28,29]. These working pairs generally have no crystal restrictions, 
and therefore, have relatively wide operating ranges. 

Hydrofluorocarbons (HFCs) are refrigerants with benign tech- 
nical properties, which are extensively used in the vapor compres- 
sion refrigeration cycle. It is worth mentioning that HFCs used 
in the absorption refrigeration cycle are almost harmless to the 
atmospheric ozone, whereas chlorofluorocarbons (CFCs) are 
harmful to the atmospheric ozone commonly cannot be used for 
the absorption refrigeration cycle [30]. HC working pairs have also 
attracted widespread attentions because of their naturally eco- 
friendly feature [8,31]. However, due to the flammability and 
explosibility, researches and applications of HC working pairs are 
impeded by numerous restrictions. 

Both H20/LiBr and NH3/H2O systems can be developed for 
future applications, but existing problems need to be solved 
urgently. Unlike H2O and NH3, absorbents suitable for alcohols, 
HFCs, and HCs are still quite few. Considering that HFCs and HCs 
can meet the refrigeration technology demands of H20 and NH3, 
absorption cycle working pairs containing HFC or HC have great 
potential for exploration and are expected to fill the gaps of H20/ 
LiBr or NH3/H20 systems. Around exploring advanced absorbents 
with HFC and HC, researchers are exerting efforts to develop novel 
absorption cycle working pairs [27,32]. 


Requests of future development 


Abate or avoid the problems of corrosion and 
crystallization 
Reduce the separating energy consumption 
1) Meet needs of the R&D of new refrigerants, 
e.g., R152a, R32 
2) Suitable for novel hybrid cycles and co- 
generation cycles 


Better techno-economic significance, and 
environmentally friendly characteristics 
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1.2. Basic properties and characteristics of ILs used as absorbent 
species 


ILs are composed of ions, and generally consist of organic 
cations containing nitrogen or phosphorus and organic or inor- 
ganic anions [33]. In 1992, Wilkes et al. [34] published reports on 
synthesis and performance of 1-ethyl-3-methylimidazolium tetra- 
fluoroborate ([EMIM]BF,). They found that this IL is moisture 
stable, does not undergo hydrolysis, has a low melting point, and 
has good thermal stability. Thereafter, the study of ILs has been 
greatly promoted, constantly developed, and emerged worldwide 
upsurge. 

Variety species of ILs exist, and their physicochemical proper- 
ties can be regulated by changing the type and structure of anions 
and cations. Thus, ILs are called designable solvents [35,36]. 
Generally, based on the types of organic cation, ILs can be divided 
into imidazolium, pyridinium, tetraalkyl ammonium, tetraalkyl 
phosphonium, and so on [37,38]. Common cationic structures 
are shown in Fig. 1. Among ILs mentioned above, imidazolium 
ILs are the most stable, most common, and most studied. 

According to the solubility differences of ILs and water, ILs 
can be divided into two categories: hydrophilic ILs, such as 
1-alkyl-3-methylimidazolium chloride ([C,MIMIJCl, n is the 
number of C atom in the alkyl, n=1, 2, 3,...), 1-alkyl-3-methyl- 
imidazolium  tetrafluoroborate ([C,MIMJBF,),  1-alkyl-3- 
methylimidazolium dimethylphosphate ([C,MIM]DMP), 1- 
alkyl-3-methylimidazolium methylsulfate ([C,MIM]MeSO,), 
and 1-alkyl-3-methylimidazolium ethylsulfate ([C,MIMIJEt 
SO,), and hydrophobic ILs, such as 1-alkyl-3-methylimidazo- 
lium hexafluorophosphate ([C,MIMJ]PF.), 1-alkyl-3-methyl- 
imidazolium trifluoromethanesulfonate ([C,MIM]TfO), and 1- 
alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide 
({C,»MIM]Tf2N) [39,40]. The hydrophilicity and hydrophobicity 
of ILs are affected by type and structure of anions and cations. 
Usually, a shorter alkyl side chain leads to stronger hydrophi- 
licity of IL, conversely, a longer alkyl side chain leads to 
stronger hydrophobicity of IL. 

As a new kind of environmentally friendly solvents, ILs have 
been widely utilized in chemical reactions, catalysis, separation, 
and electrochemistry [36,40,41]. Moreover, ILs can be used as 
absorbent species in combination with refrigerants to constitute 
new working pairs for absorption cycle. Compared with organic 
solvents, the advantages of IL are mainly embodied in the follow- 
ing aspects. 


(1) Most ILs exist as liquid at a wide range around room tem- 
perature. The boiling points of ILs are much higher than those 
of refrigerants because of the ignorable vapor pressure. In the 
generation process, the refrigerant can be easily separated 
from IL with high purity in the cycle. 

(2) The heat capacity of IL solution is small, which is beneficial in 
improving the cycle efficiency. 

(3) Solubilities with inorganic or organic species and the affinity 
with refrigerants are good, which is in favor of enhancing the 
mass transfer. 

(4) ILs have high chemical and thermal stability, high thermal 
decomposition temperature, and are non-flammable. 
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Fig. 1. Structure of the cations of the common ILs. 


(5) The chemical and physical properties of ILs can be adjusted by 
the design of anion and cation. 


Many research groups have investigated the thermophysical 
properties of working pairs consisting of ILs and common refrig- 
erants, such as H20, NH3, alcohol, and HFCs, for absorption cycle 
[42-45]. These researches exhibit the great potential and signifi- 
cant application prospects of these novel working pairs. 


1.3. Previous studies 


In 2004, Kim et al. [42] proposed some working pairs composed of 
imidazolium ILs (1-butyl-3-methylimidazolium tetrafluoroborate 
([BMIM]BF,), 1-butyl-3-methylimidazolium bromide ([BMIM]Br)) 
and refrigerants (TFE, 1,1,1,2-tetrafluoroethane (R134a)). Among these 
systems, TFE/[BMIM]Br and TFE/[BMIM]BF; were selected as alter- 
native working pairs for absorption refrigeration cycle and heat 
pumps. The investigated results of heat capacity and vapor pressure 
of the two systems show that [BMIM]Br and TFE had a stronger 
interaction, and IL working pairs generally have good thermodynamic 
properties and deserve further study. 

Since 2006, Shiflett and Yokozeki [46-48] from DuPont Com- 
pany have published lots of patents and papers on the use of ILs as 
absorbents, which are impressive researches at this field in recent 
decade. The solubilities and diffusivities of H20 [49], NH3 [23,50], 
CO, [44,51] and HFCs [52-62] in ILs were measured, and some 
new working pairs for absorption cycle were proposed, such as 
H20/[EMIM]BF,4, NH3/[DMEA]Ac (N,N-dimethylethanolammonium 
acetate), and R134a/[EMIM]BEI (1-ethyl-3-methylimidazolium bis 
(pentafluoroethylsulfonyl)imide). 

In 2006, Sen and Paolucci [45] presented that IL 1-butyl-3- 
methylimidazolium hexafluorophosphate ([BMIM]PF,) can be 
used in absorption refrigeration cycle with CO2. Some properties, 
such as vapor pressure, thermal stability, and dissolve ability, were 
discussed, which demonstrated that thermodynamic data of IL 
working pair systems are still lacking. In 2010, Martin et al. [63] 
selected appropriate IL and supercritical CO as working pairs and 
calculated the coefficient of performance (COP) for absorption 
refrigeration cycle based on the phase equilibrium of the system. 
However, because of the necessity of operating with a higher 
circulation ratio, the COP was lower than that of the conventional 
system NH3/H20. In 2012, Kim et al. [64,65] conducted theoretical 
researches about the thermodynamic performance of a miniature 
absorption refrigeration system with refrigerant/IL as working 
pairs. Previously studied ILs including  1-ethyl-3-methyli- 
midazolium __ bis(trifluoromethylsulfonyl)imide ([EMIM]Tf2N), 
[EMIM]BF4, [BMIM]BFy, [BMIM]PF.e, 1-hexyl-3-methylimidazo- 
lium _ bis(trifluoromethylsulfonyl)imide ([HMIM]Tf2N), 1-hexyl-3- 
methylimidazolium tetrafluoroborate ([HMIMJ]BF,), and 1-hexyl- 
3-methylimidazolium hexafluorophosphate ({HMIM]PF.), as well 
as refrigerants including H20, pentafluoroethane (R125), R134a, 
1,1,1-trifluoroethane (R143a), 1,1-difluoroethane (R152a), difluor- 
omethane (R32), 1,2-dichloro-1,1,2,2-tetrafluoroethane (R114), and 
1,1,2,2-tetrafluoroethane (R134). COP and feasibility were calcu- 
lated and compared. Results showed that among these systems, 
H20/[EMIM]BF, had the highest COP with a value of 0.91. 

Liang et al. [66] and Zhao et al. [67] determined the heat 
capacity, viscosity, density and phase equilibrium for methanol/ 
[BMIM]Cl (1-butyl-3-methylimidazolium chloride) and methanol/ 
[DMIM]DMP (1,3-dimethylimidazolium dimethylphosphate) sys- 
tems. They also analyzed the potential of the two systems, which 
were found to meet the requirements of benign working pairs of 
absorption cycle. 

In 2011, Zhang et al. [68] evaluated the performance of the H20/ 
[EMIM]DMP (1-ethyl-3-methylimidazolium dimethylphosphate) 
system in absorption refrigeration cycle, results indicated that its 
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COP was lower than that of the traditional H2O0/LiBr system but 
still higher than 0.7. Moreover, the generation temperature of the 
H20/[EMIM]DMP system was lower than that of the H20/LiBr 
system, showing that the new cycle could be driven by waste heat 
or hot water with lower temperature. Zhang et al. [69,70] also 
simulated the thermodynamic performance of the heat pump 
using H20/[EMIM]DMP and H20/[DMIM]DMP systems as working 
pairs, and compared the result with that of H2O/LiBr and TFE/E181 
systems. The results indicated that the COP of H20/[EMIM]DMP 
and H20/[DMIM]DMP systems were slightly lower than that of the 
H20/LiBr system, but significantly higher than the value of the TFE/ 
E181 system. 

Since 2007, author's group has conducted numerous theoretical 
and experimental studies on IL working pairs, and published 
several related papers [71-83]. Based on the solution thermody- 
namic theory, the excess Gibbs function (GĦ) evaluation method is 
expanded to IL working pairs using the UNIFAC model as calcula- 
tion tool, and this method is applied in the screening of H2O/IL, 
NH3/IL, and HFC/IL systems. According to the evaluation method, 
hydrophilic IL 1,3-dimethylimidazolium chloride ([DMIM]Cl) was 
selected and synthesized. A variety of experimental methods were 
adopted to determine thermophysical properties, such as density, 
heat capacity, vapor pressure and solubility, for H2O/IL, NH3/IL, and 
HFC/IL systems. Appropriate thermodynamic models were chosen 
to regress and predict the properties of selected systems. The 
influence of ILs as additives in traditional H20/LiBr and NH3/H20 
systems was investigated. Furthermore, the characteristics of 
absorption refrigeration cycle and heat pump using IL working 
pairs were simulated and calculated. 

Refrigerant/IL systems have increasingly received much atten- 
tion as novel, environmentally friendly working pairs. Researches 
show that IL working pairs can extend the operating temperature 
range of the cycle and strengthen the separation and generation 
processes of absorbate, thus reducing the energy consumption and 
potentially reducing or eliminating corrosion and crystallization. 
Proper combinations of refrigerant and IL may generate new 
working pairs with excellent performance. However, some pro- 
blems with IL working pairs are still encountered, such as high 
circulation ratio and low energy conversion efficiency, showing 
the necessity for further development. Fortunately, ILs have huge 
species number, and their functions can be designed by changing 
the structure of cation and the species of anion to meet the 
increasing demands for new working pairs. Therefore, selecting 
more valuable working pairs from a large number of research 
objects, measuring necessary chemical and physical properties of 
new systems with high accuracy, establishing corresponding 
models of physical and chemical properties, and evaluating the 
application significances of new IL working pairs, are still challen- 
ging issues that must be overcome to accelerate the application 
progress of IL working pairs. 


1.4. Focus and proposals of this work 


This work aims to exhibit the international research frontier of 
absorption cycle IL working pairs in recent years to make relevant 
researches more rational and efficient. The evaluation method for 
IL working pairs was introduced, which was based on macroscopic 
properties and intermolecular interactions, combining the UNIFAC 
model with extreme G? criterion. The research progress on the 
thermophysical properties of IL working pairs was presented, 
including the measurement method and modeling of vapor 
pressure, solubility, heat capacity and density for popular systems. 
Particularly, on the basis of researches on thermophysical proper- 
ties of systems containing refrigerants H20, NH3, and HFCs, several 
binary and ternary systems with potential have been presented. 
Finally, the performance of absorption cycle using IL working pairs 


was calculated and evaluated to verify further the future develop- 
mental potential of new systems. 


2. Selection of IL absorbent for working pair innovation 
2.1. General behaviors of the vapor-liquid equilibrium 


The activity coefficient of the refrigerant in IL can be calculated 
using the vapor-liquid equilibrium (VLE) data of IL working pairs 
and can be used to quantitatively characterize the volatility of 
refrigerant and the solubility in IL, establish the thermodynamic 
model of IL working pair, analyze the influence of component 
structure on the interaction between refrigerant and IL, and guide 
the selection and optimization of IL absorbent. 

For refrigerant (1)/IL (2) binary systems, the concentration of IL 
in the gas phase can be neglected, i.e., the gas phase mole fraction 
of the refrigerant species is equal to 1. Hence, the VLE relationship 
can be expressed as [84] 


p(T) = x17, P\ (T) (1) 


where p is the system pressure, x1, 71, and p§ are the liquid phase 
mole fraction, activity coefficient, and saturated vapor pressure of 
the refrigerant species, respectively. Generally, the activity coeffi- 
cient model, Wilson model, UNIQUAC model, NRTL model, and so 
on can be adopted to correlate the y4 based on the experimental 
VLE data of working pairs. Appropriate model is selected by 
considering several aspects of theoretical analysis and experimen- 
tal measurement. Therefore, the activity coefficient value of the 
refrigerant species yı can be correlated and predicted, while the 
activity coefficient value of the IL species yz cannot be calculated 
directly, because ILs have negligible vapor pressure. However, 
when y, can be estimated, y can be determined by integration 
of the Gibbs-Duhem equation at constant temperature and pres- 
sure [85]. 


xı aln y4 /0X2+X2 AlN yz /dx2 = 0 (2) 


where x is the liquid phase mole fraction of the IL species. For 
instance, the activity coefficient of refrigerant species can be 
represented by an appropriate model [85] 


In y; = 4x5 + bx3 + cx (3) 


where a, b, and c are concentration-independent parameters. The 
value of In y2 for the IL species can then be obtained in terms of 
the same parameters. 

In addition, the partial pressure of Raoult's law used for 
describing the VLE behavior of the refrigerant in a liquid ideal 
mixture is expressed as [84] 


p(T) =xip{(7) (i= 1,2,...,N) (4) 


In a p—x(y) diagram, Eq. (4) exhibits a straight line. The p—x(y) 
curve of the non-ideal mixtures that display positive deviation 
from Raoult's law falls above the line; whereas the p—x(y) curve of 
the non-ideal mixtures that show negative deviation from Raoult's 
law falls below the line, causing the negative deviation of the 
liquid phase from the ideal mixtures, exhibiting the value of 7; is 
less than 1. 

In the view of vapor pressure criterion for selecting absorption 
cycle working pairs, systems with strong absorbing ability are 
usually classified as exhibiting negative deviation from Raoult's 
law [86]. Working pairs exhibiting highly negative deviation from 
Raoult's law produce the best results, as Morrissey has described 
[87]. The solvent vapor pressure decreases because the solvent 
molecules in the liquid phase have strong affinity to the solute 
molecules. 
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Fig. 2. Vapor pressure-composition diagram of several IL working pairs. m, H20/ 
[DMIM]DMP, 323.2 K [72]; 4, TFE/[EMIM]BF,, 323.15 K [80]; e, R32/[EMIM]Tf2N, 
283.15 K [54]; ¥, NH3/[BMIM]PF¢, 283.4 K [23]. Dashed line: Raoult's law. 


Table 2 
Bond energy and bond length of some hydrogen bonds [86]. 


Hydrogen bond Bond energy/kJ mol`! Bond length/pm Compounds 


F-H...-F 28.1 255 (HF); 
O-H..--O 18.8 276 H,0 (ice) 
O-H..--O 25.9 266 CH30H, C2H50H 
N-H...--F 20.9 268 NH4F 
N-H.-.--O 20.9 286 CH3CONH; 
N—-H..--N 5.4 338 NH3 


Studies on available IL working pairs of binary and ternary 
systems and common refrigerants H20, NH3, alcohols, HCs, and 
HFCs have shown that the VLE behaviors of IL working pairs 
generally manifest as non-ideal mixtures. For example, the vapor 
pressure-composition (p—x) diagram in Fig. 2 describes the 
behavior of several IL working pairs showing negative deviation 
from Raoult's Law. 


2.2. Effects of the molecular structure 


Essentially, the non-ideal VLE behavior of the working pair 
depends on intermolecular force, which includes the interaction 
between different molecules and homologous molecules for a 
multi-species system. The quantitative relationships for contacting 
intermolecular forces and macroscopic properties are now con- 
fined to simple idealized systems. The mechanism underlying 
solubility difference among IL working pairs is not clear at the 
molecular level, the hydrogen bond is only known to play an 
important role [62,88]. 

The hydrogen bond can be simply expressed as X—H.- - - -Y, 
where the dotted line represents the hydrogen bond. Based on the 
formation theory of hydrogen bond, the strength of the hydrogen 
bond is associated with the electronegativity of X and Y atoms 
connected with the H atom. The stronger the electronegativity of 
X and Y atoms is, the stronger the hydrogen bond will be. 
Meanwhile, the hydrogen bond is also influenced by the radius 
of the Y atom. A smaller the radius of the Y atom results in a 
stronger hydrogen bond. The bond energy and bond length of 
some hydrogen bonds are listed in Table 2 [89]. Some atomic 
electronegativity values (Pauling scale) are shown in Fig. 3 [90]. In 
the H20 molecule, the electronegativity of the O atomic is larger, 
so that the H atom has excess force to form a hydrogen bond with 
highly electronegative atom in other molecules, such as F, O and N. 
Fig. 4(a) schematically shows the hydrogen bond between 


IA 
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Fig. 3. Some atoms electronegativity values (Pauling scale). 


hydrophilic IL and H20. As shown in [DMIM]CI aqueous solution, 
the H atom in H20 molecule and the Cl atom in IL anion of form 
hydrogen bond O—H . - - -Cl. Fig. 4(b) shows the hydrogen bond- 
ing interaction between R32 and [BMIM]Tf,N (1-butyl-3-methyli- 
midazolium bis(trifluoromethylsulfonyl)imide). The imidazole ring 
of IL has three hydrogen bond donors (H1, H2, H3), but only two 
proton acceptors (two F atoms) in the R32 molecule, so two 
hydrogen bonds are formed. The anion [Tf.N]~ has 10 proton 
acceptors (01, 02, 03, 04, F1, F2, F3, F4, F5, and F6), but the R32 
molecule only has two hydrogen donors (two H atoms), so two 
hydrogen bonds are also formed. 

Generally, the electronegativity of the C atomic is small, and 
C—H is generally difficult to form hydrogen bond. However, the 
C atom between two N atoms in the imidazole ring of cation will 
form a hydrogen bond because the electronegativity increases due 
to the highly electronegative neighboring two N atoms. Therefore, 
in IL aqueous solution, the O atom in a H20 molecule can easily 
form a hydrogen bond with the H atom on the C atom between 
two N atoms in the imidazole ring, while the interaction with the 
H atoms on the other two C atoms in the imidazole ring is weak. 
When the length of alkyl side chain of cation in IL decreases, the 
volume of cation decreases, the surface charge density increases, 
and the electrostatic interaction between cation and H20 molecule 
enhances. Correspondingly, the ability of IL binding together with 
H20 molecule enhances, showing that a shorter cationic alkyl side 
chain results in better affinity of IL with H20. Klimavicius et al. [91] 
studied hydrogen bond formed in D20 solutions of 1-decyl-3- 
methylimidazolium bromide ([Cj)MIMJBr) and 1-decyl-3- 
methylimidazolium chloride ({C;oMIMJCl) with 'H, °C NMR and 
Raman spectroscopy monitoring of proton/deuteron process. They 
found that anions and aggregation effects play crucial role in IL 
aqueous solutions, anions interact with cations via H-bond and 
bind water molecules in their solvation shells. 

In the study of D'Angelo et al. [92], the interaction of H20 and 
Br~ in [BMIMJ]Br was investigated combining extended X-ray 
absorption fine structure spectroscopy and molecular dynamics 
simulations. A picture of structural properties of [BMIM]Br/H20 
was gained, in the surroundings of both the imidazolium cation 
and anion Br~. This new approach has been found to be particu- 
larly well suited to obtain detailed information on IL solutions 
despite their complexity. 

For the working pairs adopting water as refrigerant, the 
interaction between water and absorbent needs to be enhanced, 
so the IL, which has good affinity with water, is studied. The 
number of ILs is beyond count, because the structure can be 
designed. The different configurations of the anion and cation 
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Fig. 4. Hydrogen bonding interaction between refrigerants and ILs. (a) H2O/[DMIM]CI and (b) R32/[BMIM]Tf,N. 


Table 3 
Relationship of IL anion and water solubility. 


Hydrophobic > Hydrophilic 
PFs BF; Brah I> 
[(CF3SO2)2N]~ CF3SO3 CH3COO~, (CH3),PO; 
[(CoF5SO2)2N]~ CH30SO; , C2H50S03 


make ILs show different hydrophilicities. The interaction between 
IL and water is influenced by the structure, including the organic 
cationic species, cationic substituent structure, and anionic spe- 
cies. Studies show that for imidazole ILs, increasing of the alkyl 
side chain of the cation leads to reduce the hydrophilicities of ILs, 
such as [DMIM]X > [EMIM]X > [BMIM]X > [HMIM]X (X means a 
certain kind of anion). 

The anion characteristics greatly affect the solubility of IL in 
water. Table 3 lists the relationship between the IL anion and 
water solubility. Generally, for the solubility of the refrigerant in IL, 
the anion shows greater effect than the cation [93], because the 
imidazolium cation only depends on the length of the alkyl side 
chain or structure to show the impact, while the anion may have 
greater variety and structural changes. 

Similar to H,O/IL systems, the hydrogen bond has a very 
important function in HFC/IL systems. The hydrogen bond heavily 
influences the physical and chemical properties of the IL working 
pairs. Forming hydrogen bond between solute and solvent will 
increase the solubility, so IL having good affinity with HFC is 
sought for HFC/IL working pairs. As in hydrophilicities, the affinity 
of IL with fluoride is also influenced by the cation species, 
structure of cation substitution, and anion species. Available 
experimental results show that for imidazolium ILs, the longer 
the cationic alkyl side chain is, the better the affinity with HFC will 
be, as shown in the following rule: [EMIM]X < [BMIM]X < [HMIM] 
X<[OMIM]X ([OMIM] means 1-octyl-3-methylimidazolium), 
which is opposite to the hydrophilic of IL [75]. 

Comparing VLE data of several HFC/IL systems, the effect of 
anion on the solubility of HFC in IL was as follows: [C,MIM]Tf2N > 
[C»MIM]PF, > [Cn»AMIM]TfO ~ [C,MIM]BF, [75]. The solubility 


increases with increasing number of fluorine atoms in the anion. 
The hydrogen bond is only one aspect of the contribution of the IL 
affinity with fluoride, far from being able to explain the pros and 
cons on the overall performance of working pairs. Thus, further 
studies are needed. 


2.3. Infinity dilution activity coefficient criterion 


The activity coefficient of the refrigerant in working pairs 
relates to the composition of the system. However, the infinite 
dilution activity coefficient of the refrigerant (y?°), calculated on 
the basis of Henry's law, is a physical quantity that is unrelated to 
the composition of the system and reflects the regular interaction 
between refrigerant and IL. At given temperature and pressure, 
Henry's constant in working pairs is given by 


kı = lim (p /X1) =y Vp; (T) ©) 


where kı is Henry's constant. At low vapor pressure of the 
refrigerant, the solubility of the refrigerant has a linear relation- 
ship with vapor pressure, and y7 and kı can be obtained when the 
concentration tends to zero, i.e., the slope of the line when 
pressure tends to zero. In fact, yf and k are direct reflections of 
the solubility of the refrigerant at low vapor pressure. Smaller 
values of yf and kı mean higher solubility of the refrigerant. 
Therefore, y7 and kı can be used to compare the degree of 
miscibility of the refrigerant and IL. 

The absorption potential was proposed as another evaluation 
criterion for absorption cycle working pairs [83]. 


w= 1/rP (6) 


where y is the absorption potential. At given temperature and 
pressure, higher values of y4 result in stronger absorbing ability of 
the absorbent. 

In Table 4, y7 decreases with increasing length of cation alkyl 
side chain. In Table 5, except for 1-butyl-3-methylimidazolium 
trifluoromethanesulfonate ([BMIMJTfO), y% decreases with 
increasing number of fluorine atoms of anion. The y? of binary 
systems composed of nine HFCs, such as R32, R134, with IL [BMIM] 
PF, are listed in Table 6. The HFC containing more carbon atoms 
and fluorine atoms shows better miscibility with [BMIM]PFg. In 
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Table 4 
Activity coefficients at infinite dilution of R32 in [PF.]-based ILs at 298.2 K. 
Item [OMIM]PF, [HMIM]PF, [BMIM]PF, [EMIM]PF, 
AO a A A AAS RAN e 
Absorbent \_J/ ne Nef rir rie VOR LF 
structure a pop poh p PISF 
F F F F 
ne 0.4525 0.5559 0.7328 1.1028 
WY 2.2099 1.7989 1.3646 0.9068 
Pi /MPa 1.6919 1.6919 1.6919 1.6919 
ky 0.7656 0.9405 1.2398 1.8658 
Table 5 
Activity coefficients at infinite dilution of R32 in the [BMIM]-based ILs at 298.2 K. 
Item [BMIM]Tf.N [BMIM]TfO [BMIM]PF. [BMIM]|BF, 
NA YA NIAS NAS 
Absorbent NG. NP O T \ ane oe 
ey H oHe S ee 
structure rl od Kp i FOISF | 
F o F F 
ne 0.5026 0.7112 0.7328 0.8794 
yi 1.9897 1.4061 1.3646 1.1371 
Pi /MPa 1.6919 1.6919 1.6919 1.6919 
ky 0.8503 1.2033 1.2398 1.4879 


Table 6 
Activity coefficients at infinite dilution of various HFCs in the IL [BMIM]PF, at 
298.2 K. 


Items R32 R134 R41 R23 R152a 
Refrigerant structure CHF, CHF,CHF, CH3F CHF3 CH3CHF, 
aie 0.7328 0.8289 0.8677 1.2847 1.3957 
ya 1.3646 1.2064 11525 0.7784 0.7165 
pi [MPa 1.6919 0.5268 3.8383 4.7044 0.5973 
kı 1.2398 0.4367 3.3305 6.0437 0.8337 
Items R161 R134a R125 R143a 
Refrigerant structure CH3CH)F CF3CH2F CHF,CF3; CH3CF3 

y 1.4969 1.8615 3.1726 3.6674 

ya 0.6680 0.5372 0.3152 0.2727 

pi [MPa 0.9230 0.6664 1.3797 1.2632 

kı 1.3817 1.2404 4.3772 4.6327 


addition, comparing R134 with R134a, y7 of R134a is much larger 
than that of R134 because of isomerization. 


2.4. VLE prediction without experimental data by the UNIFAC model 


The activity coefficient of the refrigerant species in IL working 
pairs is the basic data needed for selecting suitable IL absorbents. 
However, the obtained experimental data are sometimes insuffi- 
cient. The UNIFAC model provides a calculation method to corre- 
late and predict the activity coefficient without experimental data, 
which was proposed by Fredenslund based on the group- 
contribution method in 1975 [94]. The UNIFAC model is specifi- 
cally suitable for mixtures of molecules, but not for mixtures 
containing ionic components, because the residual term of the 
model only considers short-range interactions. To extend the 
application range of the UNIFAC model, the IL can be considered 
as a neutral molecule or a weak electrolyte solution. The reason- 
ability of this perspective can be explained from the following 
aspects [81]. The conductivity of IL is usually very low, at 
approximately 107° mS/cm. After dilution in the polar solvent 


the conductivity will increase by 10-20 times, showing strong 
association between the anion and the cation, and forming 
electroneutral ionic clusters [95]. The VLE and LLE data of IL 
systems can be well described by the GË model applied for 
traditional non-electrolyte solutions, such as the NRTL activity 
coefficient model and PR equation of state (EOS), showing that ILs 
can be reasonably regarded as molecular components and that the 
electrostatic interaction can be ignored [54,96]. Using electrospray 
ionization mass spectrometry, a large amount of ionic clusters in 
different degrees of association or supramolecular polymerization 
has been detected in ILs, showing that ILs are similar to neutral 
molecules or weak electrolytes that can be approximately treated 
as neutral components. 

In the UNIFAC model, the activity coefficient is calculated as a 
sum of two terms [97], 


Iny;= Iny$+ Iny? (7) 


where y; is the activity coefficient of component i, yF is the 
combinatorial term, and 7 is the residual term. The combinatorial 
term mainly considers the difference of molecular size and shape. 
The residual term mainly considers the interaction force between 
molecules, and can be obtained from the properties of pure 
substances and mixtures. 

To use the UNIFAC model, ILs need to be segmented into 
groups. ILs can be decomposed through the method presented 
by Kim [98] and Lei [99]. As shown in Fig. 5, the IL [BMIMICI is 
divided into one CH; group, three CH3 groups, and one [MIM]Cl 
group. Other ILs can also be divided into groups in the same way. 

Based on the study conducted by Lei et al. [99], Dong et al. [71] 
studied a series of H2O/IL binary systems, which mainly included 
three kinds of short carbon chain imidazolium cations, namely, 
[DMIM], [EMIM] and [BMIM], as well as seven kinds of anions, 
namely, Cl, Br, BF4, DMP, DBP (dibutylphosphate), EtSO4, and TfO. 
Parameters for 14 group pairs were then added to the UNIFAC 
parameter table. The proposed new group interaction parameters 
are listed in Table 7, which can be used to predict VLE and excess 
properties of systems containing H2O and imidazolium ILs at 


54 D. Zheng et al. / Renewable and Sustainable Energy Reviews 37 (2014) 47-68 


CH; 
Tn 
CH3 
h 
CH) 
PECPEDE EDP EEE EEE E PEEP ernie ia, ajn 
H3C “1, CH 
` eS Wong 


N + 


Ne? 


FILIIEIE EIT 


Fig. 5. Group segmentation of ionic liquid [BMIM]Cl. 


various compositions and temperatures. The new parameters are 
suitable for ILs based on seven kinds of anions, namely, Cl, Br, BF4, 
DMP, DBP, EtSO,, and TfO. Similarly, works have been carried out 
for HFC/IL binary systems [75] and NH3/IL binary systems [79]. The 
new group interaction parameters are also listed in Table 7, which 
can be used to extend the evaluating and predicting range of HFC/ 
IL and NH3/IL systems for developing alternative working pairs. 


2.5. Excess property criterion 


Closely associated with y;, there is the summation relationship 
of GF [84]. 


Cee x x[a(nG")/ani]p 7, 


N 
=RT > x, Iny; (8) 
i=1 i=l 


Nj 
For the refrigerant and IL binary systems, Eq. (8) can be written 
as 


GE = RT(x; In ¥1+X2 In y3) (9) 


where values of In yı and In y2 can be determined by Eqs. (1)-(3) 
using experimental data or with the help of the UNIFAC model. 

According to an analysis of the Flory-Huggins equation [101], 
when the attraction of heterogeneous molecules is stronger than 
that of homogeneous molecules, the value of Gë is less than zero. 
In other words, a negative value of GË? means the affinity between 
the species is relatively stronger. 

Thus, for a binary system at given temperature and pressure, a 
negative G? is beneficial for solute absorption, and the smaller the 
value of G is, the better the absorption effect will be. At given 
temperature and pressure, G? is a continuous function of system 
composition, which may exist more than one extreme point. The 
extreme point of G? presents the characteristic of thermodynamic 
properties of the system. Therefore, the extreme point of G} is 
presented as another selection criterion for evaluating suitable 
absorbents [102,103]. Thus, at given temperature and pressure, a 
smaller value of the extreme point of Gë for the solute-absorbent 
binary system indicates a stronger absorbing ability of the 
absorbent. 

To evaluate the solute-absorbent binary system, the criterion 
based on the analysis of the extreme point of GF is called “the 
criterion of minimum extreme value of the G'—x relationship”. 
The criterion based on the analysis of the absorption potential of 
solute is called “the criterion of maximum value of the absorption 
potential”. The former is a macroscopic trend analysis on the basis 
of relatively more data points, while the latter is a local state view 
that requires only a small amount of data. In the case pressed for 
available VLE data, the latter may be more convenient. However, 
the combination of two methods can be used to obtain more 
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Group interaction parameters for the UNIFAC model. 
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comprehensive evaluation conclusions. In addition, for the criter- 
ion of minimum extreme value of the G®—x relationship, the 
absolute value of minimum extreme value of the G? is commonly 
marked as GË ax- 

As shown in Table 8 and Fig. 6, for refrigerant R32 absorption 
systems with absorbents [BMIM]PF.e, [BMIM]BF,4, [BMIM]TfO, and 
[EMIM]Tf2N, the Ge. of R32/[EMIMJ]Tf2N binary system is at 
minimum, showing that [EMIM]Tf2N has better ability to absorb 
the refrigerant R32. 

Comparing data from the UNIFAC model with the experimental 
data listed in Table 8, the predicted trends of GË, are basically the 
same as the experimental trends of the y,/x; obtained by literature 
data. The assessment results obtained by the Gi, of the system 
agree essentially with the experimental data analysis. For R32 
species, the phase equilibrium ratio values of the R32/[EMIM]Tf2N 
system is relatively smaller, demonstrating [EMIM]Tf2N has better 


absorbing ability compared with other absorbents. 


3. Researches of H20 and IL systems 
3.1. Vapor pressure 


As presented above, excellent working pairs are generally non- 
ideal mixtures. The strong affinity of absorbent with refrigerant 
often leads to the lower saturated vapor pressure of solution than 
that of ideal solution. In other words, VLE behavior exhibits highly 
negative deviation from Raoult's law, which is used to describe the 
VLE behavior of ideal mixture. Therefore, the decrease in saturated 
vapor pressure of solution becomes a significant indicator for 
assessing the affinity between refrigerant and absorbent. 

In recent years, many researchers have investigated the VLE of 
IL working pairs systems. Table 9 summarizes the VLE measure- 
ment of H,O/IL systems in literatures. 

For the VLE measurement of H2O/IL systems, the most widely 
used method is the boiling point method. The principle is to 
precisely prepare a mixture with known composition, observe 
phase behavior in the equilibrium cell, determining the properties 
in equilibrium state, e.g., pressure and temperature. 

In 2004, Kim et al. [43] measured the vapor pressure of the 
systems H20/[BMIM]Br, H20/[BMIM]BF4, and H20/[HydeMIM]BF4 
(1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate). They 
proposed that these systems could be used in absorption heat pumps. 

Experimental data showed that the vapor pressure of water 
decreased with IL addition, but the decreasing degree differed 
according to IL type and concentration. For imidazolium based ILs, 
the interaction with H20 was influenced by the type of anion and 
the length of cationic alkyl side chain. Generally, a shorter cationic 
alkyl side chain leads to stronger hydrophilicity. Moreover, the 
influence of the anion was much more significant than that of the 
cationic side chain. 

In recent years, researchers have studied the thermophysical 
properties of H,0/phosphoric IL working pairs, such as [DMIM] 
DMP, [EMIM]DMP, 1-ethyl-3-methylimidazolium diethylphosphate 
({EMIM]DEP), 1-ethyl-3-ethyl-imidazolium diethylphosphate ({EEIM] 


Table 8 
The GE,,, for refrigerant R32 and IL binary systems at 298.20 K. 


Refrigerant IL 
T/K 
R32 [EMIM]T£,N CsH11F6N304S2 298.15 
[BMIM]PFg CsHisFeN2P 298.20 
[BMIM]BF4 CsgHisBF4N2 298.20 
[BMIM]TfO CoHi5F3N203S 298.14 


Experimental data 


DEP), and 1-butyl-3-methylimidazolium dibutylphosphate ([BMIM] 
DBP). All these ILs have strong affinity with water, especially [DMIM] 
DMP. Dong [72], He [109], Kato [112], Wang [117], Zhao [118] et al. 
measured the vapor pressure of the H20/[DMIM]DMP system at 
different temperatures, pressures, and compositions. The data pro- 
vide the foundation for calculating the absorption cycle, and con- 
firming the potential of the H20/[DMIM]DMP system to be used as 
working pairs. 

Based on the potential use of hydrophilic IL as absorbent 
species combined with HO to constitute alternative working 
pairs, Li et al. [76] chose [DMIMJ]Cl and 1-methyl-3-methyl- 
imidazolium tetrafluoroborate ([DMIM]BF,) as additives, added 
them to traditional working pairs H20/LiBr and H20O/LiCl, and 
explored the modification effect of IL on traditional working pairs. 
They determined the saturated vapor pressures of ternary systems 
at different temperatures and concentration ranges using the 
boiling point method. 

The correlation method of VLE data generally includes two 
categories. First is methods only using equations of state, such as 
PR, SRK, BWR, MH. Second is methods simultaneously using 
equations of state and the activity coefficient model, such as 
models of Margules, Wilson, NRTL, UNIQUAC, UNIFAC. 

For H2O0/IL systems, Doker and Gmehling [107] analyzed VLE 
data correlation of H2O/[EMIM]Tf2N and H20/[BMIM]Tf,N systems 
using three kinds of activity coefficient models including Wilson, 
NRTL, and UNIQUAC. Results showed that NRTL and UNIQUAC 
models, commonly used to correlate non-electrolyte solutions, are 
also suitable for IL working pairs with good accuracy. Carvalho 
et al. [120] determined VLE data of H2O combined respectively 
with 1-ethyl-3-methylimidazolium chloride ([EMIM]Cl), [BMIM]Cl, 
1-hexyl-3-methylimidazolium chloride ([HMIMJCl), and Choline 
Chloride ([N141(201)|Cl), the NRTL model was also used to correlate 
with data and a good fit was obtained. 

In addition, for binary and ternary systems containing H20 and 
IL, Wang et al. [117] and Li et al. [76] adopted the Antoine equation 
to correlate VLE data, the equation is expressed as 


4 á 
log(p/kPa)= > [A;+1000B;/(T/K —43.15)]w' (10) 
i=0 


100 


-400 
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Fig. 6. G’—xp32 for R32 systems at 298.20 K. =, [BMIMJBF,; ©, [BMIMJTfO; 4, 
[BMIM]PF,; A, [EMIM]T£,N. 


Cal. by UNIFAC 


P/MPa yili Ref. Geax /J* mol-? 
0.399 3.077 [53] —408.47 
0.400 3.367 [54] — 250.77 
0.400 3.676 [54] —93.51 
0.400 3.788 [73] — 207.90 


56 


Table 9 
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Summary of vapor pressure measurement for H2O/IL systems. 


IL Measuring method 

[BMIM]Ccl Dynamic recirculating still apparatus 
[EMIM]EtSO, Dynamic recirculating still apparatus 
[BMIM]MSO4 Dynamic recirculating still apparatus 
[EMIM]TfoN Computer-driven static apparatus 
[BMIM]Tf2N 

[DMIM]DMP Boilling point method 

[EMIM]Ac Static apparatus 

[HMIM]C1 

[DMIM]DMP Boilling point method 

[EEIM]DEP Quasi-static ebulliometer method 
[EMIM]BF4 Gravimetric method 

[BMIM]Br Boilling point method 

[BMIM]BF4 

[HydeMIM]BF,4 

[DMIM]DMP Computer-driven static apparatus 
[EMIM]Tf2N 

[BMIM]Tf2N 

[EMIM]DMP Bubble point method 

[EMIM]TfO Dynamic recirculating still apparatus 
[EMIM]TFA 

[EMIM]EtSO4 Transpiration method 

[EMIM]DMP Quasi-static ebulliometer method 
[DMIM]DMP Quasi-static ebulliometer method 
[DMIM]cl Boiling point method 

[DMIM]BF, Boiling point method 

[DMIM]DMP Quasi-static ebulliometer method 
[EMIM]DEP 

[BMIM]DBP 

[EMIM]EtSO, Boiling point method 

[EMIM]Cl Boiling point method 

[BMIM]Cl 

[HMIM]CI 

[Nineon]Cl 


[COC,MOR]|[FAP] 
[COC,MOR]|TE,N] 
[COC,PIP][FAP] 
[COC,PIP][ TEN] 
[COCPYR][FAP] 
[COC,PYR][TEN] 
[BMPIP][N(CN)2] 
[BMPYR][N(CN),] 
[BMPy][N(CN)2] 
[BMPy][SCN] 


Dynamic method 


where T is the temperature, and w is the mass fraction of the 
absorbent. A; and B; are equation parameters obtained by regres- 
sing VLE data of the system. 


3.2. Heat capacity 


Heat capacity is an important factor affecting the coefficient of 
performance for absorption refrigeration cycles. A smaller heat 
capacity benefits heat transfer, reduces energy consumption, and 
improves the cycle performance. Measuring the constant pressure 
heat capacity enables the calculation of enthalpy, entropy, Gibbs 
free energy, and other thermodynamic functions. To discuss the 
relationship between the thermodynamic model of new working 
pairs and circulation mechanism, the heat capacity of new work- 
ing pairs for absorption cycle must be studied. 

Previous studies on the heat capacity determination of IL 
working pairs in literatures mainly focused on water or alcohol 
with IL systems, whereas the heat capacity of HFC or NH3 with IL 
systems was rarely reported because of restricted test conditions. 
Table 10 summarizes researches on the heat capacity of H2O/IL 
systems published in recent years. 

Heat capacity is commonly measured by differential scanning 
calorimetry (DSC) and the micro-calorimeter method. 


Ref. Range of data 
T/K p/kPa 
104] 373.15 — 441.11 101.3 
105] 373.15 — 410.58 101.3 
106] 373.15 — 374.12 101.3 
107] 353.15 3.193 — 47.422 
353.15 5.317 — 47.156 
72] 302.57 — 441.87 2.14— 101.01 
108] 283.1 — 403.17 0.215 — 104.8 
283.15 — 423.16 0.298 — 262.8 
109] 329.35 — 429.75 6.62 —39.15 
110] 317.63 — 371.31 9.313 —57.075 
111] 373.15 — 403.15 2.5— 80.2 
43] 304.8 — 457.4 3.5— 103.2 
311.4- 475.2 3.5 — 102.6 
325.8 — 464.8 3.3 — 100.8 
1121 353.15 0.09 — 47.82 
353.15 3.19 —47.39 
353.15 5.32 — 47.16 
113] 320.35 — 411.15 6.27 — 43.56 
114] 323.3 — 363.3 5.9-69.3 
328.3 — 368.2 6.3 — 84.4 
115] 302.9 — 322.9 0.345 — 12.141 
116] 323.76 — 367.26 12.609 — 58.804 
117] 328.936 — 379.411 13.834 — 101.270 
80] 287.15 — 437.45 1.10 — 100.98 
82] 312.25 — 403.60 4.69— 103.23 
118] 312.88 — 373.01 6.937 — 79.824 
324.01 — 373.26 11.872 — 83.919 
317.98 — 368.07 9.184- 81.630 
119] 320.35 — 373.35 4.21 — 24.75 
120] 355.58 — 435.03 50—100 
354.92 — 430.10 
354.64 — 405.30 
355.11 — 411.12 
121] 290.1 — 351.7 100 
291.0 — 349.6 
301.2 — 349.3 
313.2 — 353.9 
307.4- 351.6 
291.8 -357.5 
237.64— 273.15 
236.26 — 273.15 
244.19 — 300.30 
238.03 — 273.15 
Table 10 


Summary of heat capacity measurement for H2O/IL systems. 


System Measuring method Ref. Range of T/K 

H20/[DMIM]DMP Calvet heat flow calorimeter 72)  303.15—353.15 
H20/[EMIM]EtSO, DSC 122] 283.15—343.15 
H20/[EMIM]TfO 283.15 — 343.15 
H20/[EMIM]TFA 283.15 — 343.15 
H20/[EMIM]EtSO,4 DSC 123] 293.15—318.15 
H20/[BMIM]CH3SO4 293.15 — 318.15 
H20/[EMIM]TfO 293.15 — 318.15 
H,0/[BMIM]TfO 293.15 — 318.15 
H,0/[DMIM]DMP Adiabatic solution calorimeter [109] 298.15—323.15 


H20/[EMIM]EtSO4 DSC 124 303.2 — 353.2 
H20/[EMIM]TfO 303.2 — 353.2 
H20/[BMIM]CH3SO, DSC 125) 303.2 — 353.2 
H,0/[BMIM]TfO 303.2 — 353.2 
H20/[BMIM]BF4 DSC 126] 278.15—333.15 


H30/[EMIM]DMP 
H30/[EMIMJEtSO, 


Adiabatic solution calorimeter [113] 298.15—323.15 
Isoperibol solution calorimeter [119] 303.15—323.15 


Heat capacity data in literatures are mostly correlated with 
temperature by a polynomial equation. A polynomial equation is 
then used to correlate the relationship of heat capacity with 
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temperature and composition. 
3 : 
Cp = > (a; +b;T)x' (11) 
i=0 


where C, is the heat capacity. a; and b; are equation parameters 
obtained by regressing experimental data. 


3.3. Density 


As an important part of the fluid PVT data, liquid density plays a 
significant role in thermodynamics and is necessary for calculating 
other physical properties and studying of heat and mass transfer. In 
practice, the density of working pairs is closely related to temperature 
and composition, so the influence of the density change of a solution 
cannot be ignored when analyzing the absorption refrigeration cycle. 

The density of most imidazolium based ILs roughly ranges 
within 1.1 — 1.6 g/cm?, which is more important than that of water 
in practical applications. The density of IL mainly depends on the 
type of anion and cation, and the influence of anion is more 
obvious. Generally, a larger anion size leads to greater IL density, 
and a larger volume of organic cation leads to smaller IL density. 

Previously studies on the density determination of IL working pairs 
in literatures mainly focused on water or alcohol with IL systems, 
whereas the density of HFC or NH3 with IL systems was rarely 
reported because of restricted test conditions. Researches on the 
density of H2O/IL systems in recent years are summarized in Table 11. 

Density is commonly measured by the pycnometer [127,128], 
gravity balance (Westphal) [109,119,129,130], and vibrating 
U-tube densimeter [74, 131—139]. The vibrating U-tube 


Table 11 

Summary of density measurement for H2O/IL systems. 
System Measuring method Ref. Range of T/K 

H,0/[BMIM]N(CN)2 

H20/[BMIM]C(CN)3 


Vibrating tube densimeter [131] 278.15 — 363.15 


278.15 — 363.15 


H20/[DMIM]Cl Vibrating tube densimeter [74] 293.15 — 318.15 
H,0/[BMIM]Br Vibrating tube densimeter [132] 298.15 
H20/[DMIM]Br 298.15 


H20/[EMIMJEtsO, 
H30/[BMIM]CH;SO, 
H30/[EMIM]TfO 
H20/[BMIM|TfO 


293.15 — 318.15 
293.15 — 318.15 
293.15 — 318.15 
293.15 — 318.15 


Vibrating tube densimeter [123] 


H,0/[BMIM]Cl Vibrating tube densimeter [133] 298.15 — 343.15 
H0/[HMIM]Cl 298.15 — 343.15 
H20/[OMIM]Cl 298.15 — 343.15 


H20/[EMIMJEtSO, 
H30/[BMIM]CH;SO, 
H,0/[DMIM]DMP 


298.15 — 328.15 
298.15 — 328.15 
298.15 — 323.15 


Vibrating tube densimeter [134] 
Vibrating tube densimeter [135] 
Gravity balance [109] 


H20/[BMIM]Cl Density bottle [136] 298.15 
H0/[EMIM]Br 298.15 
H,0/[BMIM]Br 298.15 
H20/[BMIM]BF4 298.15 
H,0/[EMIMJEtSO, Westphal balance [129] 278.15 — 333.15 


H,0/[EMIM]BF,4 Vibrating tube densimeter [137] 298.15 
H20/[BMIM]BF4 298.15 
H20/[HMIM]BF4 298.15 


H30/[EMIMJEtSO, 
H20/[EMIM]TfO 
H20/[EMIM]TFA 


278.15 — 348.15 
278.15 —348.15 
278.15 — 348.15 


Vibrating tube densimeter [138] 


H,0/[EMIM]BF, Vibrating tube densimeter [95] 298.15 
H0/[BMIM]BF4 298.15 
H20/[HMIM]BF4 298.15 
H20/[BMIM]PFg 298.15 
H0/[BMIM]DCA 298.15 


H20/[EMIM]BF4 
H20/[EMIMJEtSO, 
H30/[EMIMJEtsO, 


298.15 
298.15 
[130] 278.2-338.2 


Vibrating tube densimeter [139] 


Westphal balance 


H20/[EMIM]BF4 Pycnometer [127] 293.15 — 323.15 
H20/[BMIM]BF4 Pycnometer [128] 303.15 — 353.15 
H20/[EMIM]EtSO4 Gravity balance [119]  303.15-323.15 


densimeter method is one of the most common and highest 
accuracy methods to date. The uncertainty can reach 
+1x10~°gcm~?, and the measurement error does not exceed 
0.005%. The working principle is to change the inherent oscillation 
frequency of the vibrating tube with the mass change (i.e., the 
density change) of the liquid flowing through the vibrating tube. 
After determining the relation between the oscillation period of 
the vibrating tube and liquid density, the density can be directly 
measured. 

Solution density can be correlated with composition and 
temperature by various expressions. Given that the density of 
the IL solution has a simple relation with temperature and 
composition, a polynomial equation is selected to regress the 
density data of H2O/IL systems as 


i i i 
p= X Aw"+T Y Byw"+T? X Caw" (12) 
n=0 n=0 n=0 


where p is the density of the binary system, w is the mass fraction 
of water. An, Bn and Cn are polynomial parameters obtained by 
fitting experiment data through the least square method. 


3.4. Other properties 


For a detailed investigation of the suggested working pairs, 
more properties such as viscosity, surface tension, and thermal 
conductivity are required. 

The viscosity of IL working pair systems was rarely reported in 
literatures. The studied systems are mainly H20/[EMIM]BF,, H20/ 
[BMIM]BF4, H20/[HMIM]BF4, H20/[BMIM]Cl, H20/[EMIM]Br (1-ethyl- 
3-methylimidazolium bromide), H20/[BMIMJ]Br, H20/[EMIM]EtSO, 
(1-ethyl-3-methylimidazolium ethylsulfate), H20/[BMIM]MeSO, (1- 
butyl-3-methylimidazolium methylsulfate), and so on. Results show 
that the viscosity of H2O/IL systems gradually decreases with the 
increasing water content. The viscosity of IL is generally tens to 
hundreds times higher than that of water mainly because of the 
different degrees of IL aggregation, due to the strong hydrogen bond 
between anion and cation. When water is added to IL, the hydrogen 
bond between the anion and cation weakens, which increases the 
mobility of ions and decreases the viscosity of IL. 

The surface tension of IL is between that of water and common 
solvent, and decreases with the increase of temperature. However, 
the surface tensions of the refrigerant/IL systems are still rarely 
reported. The systems conducting research on surface tension 
include H20/[EMIM]BF4, H20/[BMIM]BF,, H20/[HMIM]BF,4, H20/ 
[BMIM]Cl, H20/[BMIM]Br, H20/[EMIM]Br, and so on. The addition 
of IL can reduce the surface tension of the refrigerant quickly, 
which is beneficial to application in absorption refrigeration cycle. 

The reports on the thermal conductivity of IL working pair systems 
which is necessary for the design of heat-transfer equipment are 
limited. The thermal conductivities for the binary systems H20/ 
[EMIMI]EtSO,4, H20/[DMIM]DMP, CH30H/[DMIM]DMP, and C2H50OH/ 
[EMIMIJEtSO, were reported. It is expected that refrigerant/IL systems 
could be perfect heat-transfer fluids. 

Moreover, properties such as toxicity, environmental impact, and 
cost have to be evaluated and considered for further investigation. 


4. Researches of other refrigerant and IL working pair systems 
4.1. NH; and IL systems (solubility and other properties) 


Because of the good mutual solubility of NH3 and H20 and the 
great latent heat of evaporation, the NH3/H2O working pair is 
widely used in various kinds of absorption refrigeration cycles. 
However, the system has some drawbacks such as high working 
pressure, toxicity, and difficulty in regenerating NH3 from H20. 
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Table 12 
Summary of solubility measurement of NH3/IL systems. 
Ionic liquid Measuring method Ref. 
[BMIM]PFe 
[HMIM]CI 
[EMIM]Tf>N 
[BMIM]BF,; 
[EMIM]Ac 
[EMIMJEtSO, 
[EMIM]SCN 
[DMEA]Ac 
[EMIM]BF4 
[BMIM]BF, 
[HMIM]BF4 
[OMIM]BF4 
[DMIM]DMP 
[BMIM]Zn2Cls 


Static method [23] 


Static method [50] 


Static method [140] 


Isothermal synthetic [78] 
Static method [141] 


Therefore, a great number of new absorbents including IL have 
been studied as alternatives. 

To find novel working pairs containing NH3 as refrigerant and IL as 
absorbent, researchers have carried out related works. The research 
progress of VLE of NH3/IL systems is presented in Table 12. The 
common method of determining the VLE data of NH3/IL systems is as 
follows. First a certain mass of IL is added to the static equilibrium cell 
that is then filled with a certain mass of NH3. Second the equilibrium 
cell is placed in a thermoregulated system (oil or water bath). The 
temperature and pressure are recorded when the NH3/IL system 
reaches equilibrium at constant temperature. Third the equilibrium 
pressures at different temperatures are recorded by varying the 
temperature. Finally an appropriate model is used to calculate and 
obtain p—T—x data. Using these experimental data, researchers can 
then verify the applicability of the thermodynamic model to NH3/IL 
working pairs, simulate the working condition of the absorption 
cycle, and explore the prospects of the working pairs in practical 
applications. 

As shown in Table 12, Yokozeki et al. [23,50] pioneered works 
in this region. They found that solubility of NH3 in ILs such as 
[HMIM]CI and [BMIM]PF, at 298 K is almost equivalent to that in 
water, i.e., negative deviations from Raoult's law for the two 
systems are at relatively the same level. Another important result 
is that NH3 does not react with ILs under the experimental 
condition, and no loss of solvent occurs even in the acid solvent 
[HMIM]CI. However, no absorption mechanism between NH; and 
various ILs is proposed and further study is deserved at the 
molecular level by spectroscopy methods. 

In another work, Yokozeki [50] calculated residual properties 
such as residual enthalpy, residual entropy, and residual Gibbs 
energy of NH3/IL systems. They found that these residual proper- 
ties are all negative. The negative value of the NH3/[DMEA]Ac 
system was the maximum, and residual properties of this system 
were close to those of the NH3/H20 system, indicates that a strong 
interaction like hydrogen bond exists between NH; and ILs. 

In 2010, Li et al. [140] determined the VLE data of NH3 and four 
ILs [C,MIM]BF, (n=2, 4, 6, and 8) and studied the effect of alkyl 
side chain length of ILs on NH3 absorption. Results show that the 
solubility of NH3 in studied ILs is high and increases with the 
increase of alkyl side chain length. They believed that increase of 
carbon chain length leads to the decrease of ILs density, generating 
more free volume, and dissolving more NH3 molecules. Using the 
Krichevisky-Kasarnovsky equation, Li et al. [140] calculated the 
Henry's constant and partial molar volume of NH3, they estab- 
lished a thermodynamic model to calculate the enthalpy, entropy, 
Gibbs energy, and heat capacity of the system. Results confirm the 
non-spontaneity of the absorption process and the disadvantage of 
high temperature in NH3 absorption. 


Range of T/K Range of p/MPa Range of x/mol% 


283.4—355.8 0.138 — 2.700 23.9 —86.2 
283.1 — 347.9 0.044 — 2.490 6.0—83.7 
283.3 — 347.6 0.114- 2.860 4.5 — 94.8 
282.2 -355.1 0.091 — 2.570 6.8 — 84.4 
282.5 —348.5 0.321 — 2.891 47.3 — 87.7 
282.7 — 372.3 0.287 — 4.777 42.4- 87.5 
283.2 — 372.8 0.244- 5.007 34.0 — 87.6 
283.2 — 372.8 0.136 — 4.249 45.4- 86.5 
293.15 — 333.15 0.11 — 0.63 11.85 — 69.21 
293.15 —333.15 0.07 — 0.83 6.08 — 75.31 
293.15 — 333.15 0.14—0.71 12.80 — 75.43 
293.15 — 333.15 0.10—0.61 13.21 — 80.81 
293.15 — 333.15 0.053 — 0.665 7.0 —75.2 
323.15 —563.15 0.067 — 1.965 83.62 —94.67 


Chen et al. [141] measured total pressures of NH3 and metal 
ion-containing IL [BMIM]Zn.Cl, by a static method, and correlated 
experimental data by the modified UNIFAC model. The solubilities 
of NH; in [BMIM]Zn,Cl, is higher than that in the other ILs tested 
by Yokozeki et al. [50], but lower than that in ZnCl. Author's group 
[78] chose NH; and [DMIM]DMP as working pairs and determined 
solubilities by the isothermal synthesis method. The solubility 
decreases with increased temperature and increases with 
increased pressure. The experimental data are consistent with 
the calculated results from the NRTL model, with a relative 
deviation within 5%. 

The literature survey reveals that properties of NH3/IL working 
pairs have been gradually examined in recent years. Given the 
controllability of the IL structure (i.e., various ILs can be synthe- 
sized according to a given target or practical needs), the gap in this 
applied research field urgently needs to be filled. Moreover, 
considering the high solubility and low energy consumption in 
the separation of NH3/IL systems, it is necessary to study the 
interaction between NH3 and various ILs, for the development of 
suitable working pairs for absorption cycle. 


4.2. HFCs and IL systems (solubility and other properties) 


As a kind of refrigerant, HFCs are also concerned in developing 
novel absorption IL working pairs. The research progress of the 
VLE of HFC/IL systems is summarized in Table 13. The studied 
refrigerants are tetrafluoromethane (R14), trifluoromethane (R23), 
R32, fluoromethane (R41), R125, R134, R134a, R152a, and fluor- 
oethane (R161). 

Shiflett and Yokozeki et al. [47,53,54,58,61,62] determined the 
solubility and diffusivity of HFCs such as R32, R134a, R125, R23, 
R152a, R134, and R161 in ILs by gravimetric microbalance, and 
used the EOS, NRTL, and Antoine-type equation to process the 
experimental data. Ren et al. [96] determined the solubility of 
R134a in [HMIM]-based ILs, regressed the experimental data using 
the PR EOS with van der Waals two-parameter mixing rules, and 
studied the effect of the anion and cationic carbon chain length on 
the solubility and molar volume of R134a in ILs. Author's group 
[73] determined the solubility of R32 and R152a in 1-ethyl-3- 
methylimidazolium trifluoromethanesulfonate ([EMIMJ]TfO) and 
[BMIM]TfO by the isothermal synthesis method, results well fit 
the NRTL model. 


4.3. HC and IL systems (solubility and other properties) 
Many researchers have studied the physical and chemical 


properties of HC and IL binary systems to discuss their potential 
use as working pairs of absorption cycle. Lee et al. [144] measured 
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Table 13 
Summary of solubility measurement of HFC/IL systems. 
Refrigerant Absorbent Measuring method 
R32 [EMIM]TfO 
[BMIM]TfO 
R152a [EMIM]TfO 
[BMIM]TfO 
R14 [HMIM]Tf,N 
R134a [HMIM]Tf2N 
[HMIM]PF, 
[HMIM]PF, 
R23 [EMIM]PFg 
R32 [BMIM]PFe 
[BMIM]BF, 
R125 [BMIM]PFg 
R134a [BMIM]PFe 
R23 [BMIM]PF, 
R152a [BMIM]PF, 
R32 [DMPIM]TMeM 
[EMIM]BEI 
[DMPIM]BMel 
[EMIM]BMel 
[Pmpy]BMel 
[bmpy]BMel 
[BMIM]Ac 
[BMIM]SCN 
[BMIM]MeSO, 
[EMIM]TFES 
[BMIM]TFES 
[HMIM]TFES 
[DMIM]TFES 
[BMIM]HFPS 
[BMIM]FS 
[BMIM]TPES 
[BMIM]TTES 
R134a [EMIM]BEI 
[BMIM]HFPS 
[BMIM]TPES 
[BMIM]TTES 
[6,6,6,14-P]TPES 
[4,4,4,14-P]HFPS 
R41 [BMIM]PFe 
R161 [BMIM]PFe 
R134 [BMIM]PFe 
R134 [EMIM]Tf.N 
R125 [EMIM]Tf2N 
R134a [EMIM]Tf2N 


Isothermal synthetic 


Synthetic method 
Static method 


Synthetic method 


solubility data of various HCs in [BMIM]Tf2N by the saturation 
method and obtained Henry's constant through the experimental 
data. Shiflett et al. [47,60,61] determined the solubility data of a 
series of HCs and fluorinated benzene in [EMIM]Tf2N using the 
gravimetric microbalance and volumetric method, then experi- 
ment data were correlated with the NRTL model. Liu et al. 
[145,146] predicted solubilities of small HCs in about 700 ILs, 
and reported solubilities of methane, ethane, ethylene and pro- 
pane in trimethyloctylphosphonium _bis(2,4,4-trimethylpenty]) 
phosphinate ([P8111]TMPP), tetrabutylphosphonium bis(2,4,4-tri- 
methylpentyl)phosphinate ([P4444]TMPP) and  [EMIMJTf2N. 
Table 14 summarizes studies on the VLE of HC/IL systems pub- 
lished in recent years. 


4.4, Alcohol and IL systems (solubility and other properties) 


Kim et al. [42] determined the solubility data of TFE in 
[BMIM]Br and [BMIMI]BF, by the boiling point method and 
correlated the data with the Antoine-type equation. Shiflett 
et al. [56] determined the VLE data of the n-butanol/[BMIM]PF. 
system using the volumetric method. Verevkin [148], Zhao 
[118], Wang [116], Jiang [110], Wang [80,117], Shen [149] and 
Carvalho et al. [120] determined the VLE data of alcohols and 
benzene in different ILs by the static method, correlated the 


Gravimetric microbalance 


Gravimetric microbalance 


Gravimetric microbalance 


Gravimetric microbalance 


Gravimetric microbalance 
Gravimetric microbalance 


Ref. Range of T /K Range of p/MPa 
[73] 273.14— 298.15 0.097 — 0.857 
273.13 — 298.19 0.104— 0.902 
273.16 — 348.17 0.040 —0.848 
273.22 — 348.15 0.056 —0.879 
[142] 293.3 — 413.3 1.202 — 9.582 
[96] 398.15 — 348.15 0.042 — 4.320 
0.081 — 14.330 
0.081 — 12.880 
[143] 308.17 — 367.33 1.602 —51.64 
[54] 283.2 — 348.2 0.0097 — 0.9999 
283.0 — 348.2 0.0097 — 0.9999 
283.1 — 348.3 0.0099 — 0.9998 
283.0 — 348.2 0.0097 — 0.3500 
282.6 — 348.1 0.0098 — 2.0002 
283.1 — 348.2 0.0097 — 0.4505 
[53] 283.15 — 348.15 0.0094 — 1.0005 
283.15 — 348.05 0.0096 — 1.0005 
298.15 0.0099 — 1.0011 


283.15 — 348.05 
283.15 — 348.05 
298.15 


0.0096 — 1.0005 
0.0095 — 1.0004 
0.0096 — 1.0000 
0.0099 — 1.0004 
0.0095 — 0.9992 
0.0099 — 1.0006 
0.0099 — 1.0016 
0.0097 — 0.9989 
0.0099 — 0.9980 
0.0096 — 1.0010 
0.0095 — 1.0004 
0.0100 — 1.0005 
0.0095 — 0.9994 
0.0095 — 0.9992 
0.0103 — 0.3505 
0.0099 — 0.3506 
0.0102 —0.3505 
0.0102 —0.3505 
0.0098 — 0.3504 
0.0099 — 0.3504 
0.0097 — 1.9995 
0.0099 — 0.7005 
0.0101 —0.3505 


[62] 283.10 — 348.10 
283.10 — 348.10 
283.05 — 348.10 
283.10 —348.15 
282.90 — 348.10 
283.05 — 348.10 
[58] 283.09 — 348.18 
283.06 — 348.16 
283.11 — 348.16 


[47] 282.9 — 348.1 0.0099 — 0.3505 
[61] 283.1 — 348.2 0.0100 — 0.9998 
282.7 — 348.1 0.0100 — 0.3503 


data using NRTL, Antoine-type and other equations, and 
obtained the activity coefficients and model parameters of 
different solutes in ILs. González et al. [150] determined vapor 
pressure of four alcohols in 1-hexyl-3-methylimidazolium tri- 
fluoromethanesulfonate ([HMIM]TfO) using the vapor pressure 
osmometry technique and derived other thermodynamic prop- 
erties. Particularly, Zhang et al. [151] proposed a new alter- 
native working pairs for absorption refrigeration cycle, which 
was ternary solutions H20/CH30OH(or C,H;OH)/[BMIM]DBP, and 
measured vapor pressure of binary solutions H20(CH30H, 
orC;Hs5OH)/[BMIM]DBP with an inclined boiling apparatus to 
obtain binary interaction parameters of NRTL activity coeffi- 
cient model, then measured and predicted vapor pressure of 
proposed solutions with the same apparatus and model. Studies 
on the VLE, density, heat capacity, and viscosity of alcohol and 
IL systems published in recent years are listed in Tables 15-18. 


5. Assessment of absorption cycle adopting IL working pairs 


5.1. Simulation for single-effect absorption cooling cycle 


To determine whether a working pair has the potential for 
further research and development, it should be placed in an 


60 


Table 14 


D. Zheng et al. / Renewable and Sustainable Energy Reviews 37 (2014) 47-68 


Summary of solubility measurement of HC/IL systems. 


Refrigerant Absorbent Measuring method Ref. Range of T /K Range of p /kPa 
Propane [BMIM]Tf2N Saturation method [144] 279.98 — 339.97 88.3 — 1062.2 
Propene 999 — 1218.2 
Butane 280.00— 340.00 334-3041 
1-Butane 21.0—3511 
Hexane [BMIM]TfO Static apparatus [147] 363.15 15.79 — 165.7 
1-Hexane 22.54 — 180.9 
2,2,4-Trimethylpentane 11.86 — 65.21 
1-Nonene 2.599 — 16.34 
Decane 2.000 —6.202 
1-Decane 1.132 — 6.989 
1-Hexane [OMIM]TfO 16.19 — 198.8 
1-Octene 1.672 — 30.24 
1-Nonene 0.784 — 13.45 
Fluorinated benzenes [EMIM]Tf2N Volumetric method [60] 282.5 — 373.0 — 
Benzene [BMIM]Tf2N Static method [148] 298.15 — 313.15 0-24.331 
Methane [P8111 ]TMPP Isochoric saturation method [145,146] 299 — 323 1090 — 3610 
Ethane 540 — 1690 
Ethylene 1080—3540 
Propane 50—210 
Methane [P4444]TMPP 313 -353 0-5000 
Ethane 
Ethylene 
Propane 
Methane [EMIM]Tf2N 299-354 0— 4000 
Ethane 
Ethylene 
Propane 
Table 15 
Summary of vapor pressure measurement of alcohol/IL systems. 
Refrigerant Absorbent Measuring method Ref. Range of T /K Range of p/kPa 
Ethanol EMIM]TOS Ebulliometric method [152 373.15 0.0533 —0.2250 
Methanol EEIM]DEP Quasi-static method [110 299.20 — 366.19 14.336 — 65.452 
Ethanol EEIM]DEP 309.80 — 370.80 9.139 — 56.716 
Methanol EMIM]DMP Bubble point method [113 301.95 — 386.15 0.01134 — 0.04630 
Ethanol 309.05 — 389.15 0.00747 — 0.05110 
Methanol DMIM]Ccl Quasi-static method [149 298.85 — 348.63 16.65 — 102.57 
Ethanol DMIM]CI 312.19 — 361.04 11.54— 101.90 
1-Butanol BMIM]PF, Volumetric method [56] 285.7 — 332.0 - 
Methanol BMIM]Tf2N Static method [148 298.15 — 313.15 0—35.450 
Ethanol BMIM]Tf2N 0— 17.928 
Propanol BMIM]Tf2N 0— 6.986 
Methanol EMIM]DMP Quasi-static method [116 299.20 — 343.86 15.969 — 60.419 
Ethanol EMIM]DMP 309.95 — 350.67 14.165 — 60.073 
1-Propanol DMIM]DMP Quasi-static method [117 330.478 — 401.916 11.094 — 100.935 
2-Propanol DMIM]DMP 316.402 — 384.177 11.149 — 102.082 
TFE BMIM]Br Boiling-point method [42] 293.2 — 469.4 0.0049 — 0.1010 
BMIM]BF,4 315.2 — 469.4 0.0068 — 0.1003 
EMIM]BF4 Boiling-point method [80] 283.15 — 460.85 0.00306 — 0.10102 
Methanol DMIM]DMP Quasi-static method [118] 285.96 — 343.71 7.811 — 76.224 
EMIM]DEP 285.31 — 342.19 7.943 — 79.563 
BMIM]DBP 286.82 — 335.10 8.431 — 80.194 
Ethanol DMIM]DMP 298.62 — 359.57 7.853 — 74.249 
EMIM]DEP 302.30 — 356.00 9.068 — 81.091 
BMIM]DBP 298.42 — 355.12 7.606 — 77.658 
Ethanol EMIM]Ccl Boilling point method [120] 334.80 — 392.42 50—100 
BMIM]CI 335.03 — 405.89 
HMIM]Cl 335.15 — 411.86 
NinonyICl 334.99 — 353.06 
H20/CH30H BMIM]DBP Boiling-point method [151] 315.95 — 392.65 11.34—72.02 
H20/C2H50H 321.65 — 395.65 11.76 — 69.02 
1-Propanol HMIM]TfO Vapor pressure osmometry technique [150] 323.15 10.39 —12.07 
2-Propanol 20.12 —23.52 
1-Butanol 3.89— 4.46 
2-Butanol 8.85 — 10.56 
1-Pentanol 1.44-1.63 


absorption cycle to assess the performance by simulation or 
experiment. Some researchers have carried out such works, but 
the main method is still limited to simulation and analysis on 


the basis of thermophysical properties of the IL working 
pairs. The works mainly focus on single-effect absorption 
cooling cycle. 
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Table 16 
Summary of density measurement of alcohol/IL systems. 


Refrigerant Absorbent Measuring method Ref. Range of T/K 
Ethanol DMIM]JDMP Gravity balance 109 298.15 — 323.15 
Methanol 
Methanol OMIM]Ccl Vibrating tube densimeter 153 298.15 — 328.15 
Methanol EMIM]DMP Vibrating tube densimeter 154 293.15 — 333.15 
BMIM]DMP 
Ethanol EMIM]DMP 
BMIM]DMP 
Ethanol EMIM]HS Viscodensimeter 155 288 — 318 
EMIM]OS 
EMIMJES 
EMIM]BS 
Ethanol EMIM]Tf£,N Density meter 156 293.15 — 323.15 
Table 17 
Summary of heat capacity measurement of alcohol/IL systems. 
Refrigerant Absorbent Measuring method Ref. Range of T/K 
Ethanol [DMIM]DMP Adiabatic solution calorimeter [109] 298.15 — 323.15 
Methanol 
Ethanol [EMIM]DMP Adiabatic solution calorimeter [113] 298.15 — 323.15 
Methanol 
Methanol [HMIM]BF4 DSC [157] 283.15 — 323.15 
Methanol [OMIM]BF4 
Table 18 
Summary of viscosity measurement of alcohol/IL systems. 
Refrigerant Absorbent Measuring method Ref. Range of T/K 
Ethanol y ; 
Methanol [DMIM]DMP Capillary viscometer [109] 298.15 — 323.15 
Methanol [OMIM]Cl Capillary microviscosimeter [153] 298.15 — 328.15 
Methanol pl pai 
ae [EMIM|DMP Micro viscometer [154] 293.15 — 333.15 
[BMIM]DMP 
Ethanol [EMIM]Tf,N Viscosity meter [156] 293.15 — 323.15 


m, [BMIM]BF,; o, [BMIM]TfO; 4, [BMIM]PFs; A, [EMIM]T£,N 


r O 
Pe 


Qe 


5 


4 


Fig. 7. Schematic of the single-effect absorption cooling cycle. 


The single-effect absorption cooling cycle consists of an 
evaporator E, a condenser C, a generator G, an absorber A, and a 
solution heat exchanger SHX, whose schematic is shown in Fig. 7. 
1—2-—5-—3-1 represents the refrigerant circulation loop, and 
6—2—5-—4-—6 is the solution circulation loop. In the cycle 
calculation, mass transfer, heat transfer, and resistance loss in 
the solution flow process are neglected. Therefore, both positive 
and reverse cycles can be considered as ideal cycles, with the 
highest efficiency of heat transfer and cooling properties. 


In the refrigerant circulation loop, the generator is heated by Qg 
and separates high-pressure refrigerant vapor 5. The vapor refrig- 
erant then enters the condenser, liquefies into liquid refrigerant 3, 
and releases heat Qc. The high-pressure liquid refrigerant passes 
through an expansion valve and goes into the evaporator. In the 
evaporator, the liquid refrigerant is vaporized and the cooling load 
Qer is produced. 

In the solution circulation loop, the low-pressure vapor refrig- 
erant 1 from the evaporator is absorbed in the absorber by the rich 
solution 8 concentrated in the generator and then releases heat Qa. 
The pressure of the poor solution 2 leaving the absorber is 
enhanced after pumping, and then the poor solution is supplied 
to the generator. In the generator, the poor solution is heated up to 
boiling by the low-grade heat source, and the refrigerant is 
vaporized and separated from the solution. The rich solution 
4 desorbing the refrigerant reduces the pressure through the 
expansion valve and flows back to the absorber. 

To simplify the computer simulation of a cycle, following 
assumptions are made. The pressures of the condenser and the 
generator are equal (pc=pc), and the pressures of the absorber 
and evaporator are equal (pa=pr). The connecting pipeline has no 
pressure drop. Compared with other forms of thermal power, the 
pumping work is so small that is usually ignored. The refrigerant 
expansion process from the condenser to the evaporator is 
isenthalpic, and the temperature of evaporator outlet is the dew 
point of the pure refrigerant. The flow rate of vapor refrigerant is 
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set to 1kgs7! 


[158,159]. 

According to the above assumptions, if the refrigerant vapor 
generated in the evaporator is 1 kgs~—'!, the heat added to the 
evaporator is (h3+ Qpr), and the heat leaving the evaporator is hp. 
Based on the energy balance, the cooling capacity Qeg of single- 
effect absorption cooling cycle is 


Qe =hz—hy (13) 


where h; is the enthalpy of stream 1 in Fig. 7, hz is the enthalpy of 
the saturated refrigerant vapor at evaporation pressure. 

In the generator, when per unit mass (e.g., 1 kg) of refrigerant 
vapor is produced, fkg poor solution is needed to flow into the 
generator. The value of fis called the circulation ratio. Conversely, 
if the poor solution flowing into the generator is fkg, 1kg 
refrigerant vapor is produced in the generator. Therefore, the rich 
solution leaving the generator is (f—1) kg. The pure absorbent 
almost has no vapor pressure, i.e., the pure absorbent in the 
stream inflow and outflow of the evaporator remains the same 
based on the mass balance, 


, and the vapor of absorbent is ignored, etc. 


fwp = (f- Dwr (14) 
or 
f = Wr/(Wr — Wp) (15) 


where wy, is the mass fraction of the absorbent in poor solution, 
w, is the mass fraction of the absorbent in rich solution. 

Investigating the heat into and out of the generator, the heat 
added to the generator is Q¢+/fh3, and the heat leaving the 
generator is (f—1)h4+hs, thus 


Qs =(f— ha +hs —fh3 (16) 


The performance of absorption refrigeration cycle is character- 
ized by the coefficient of performance defined as 


COP = Qe/Qg = (h2 — hy)/fhs +f (h4 — h3) — h4] (17) 


where h; represents the enthalpy of stream i (kJ kg~'). The 
enthalpy of each state point can be calculated on the basis of 
thermodynamic properties of refrigerant or IL working pair. 

The above analysis reveals that COP determination is based on 
the enthalpy of stream. In principle, the enthalpy of stream is the 
state function on the temperature, pressure, and composition of a 
given stream. In fact, Fig. 7 is also a phase diagram of the binary 
system. Using Fig. 7, the calculation method can be explained by 
the fact that the state of each point in the cycle needs to be 
determined, and then the enthalpy can be calculated based on the 
state. For example, based on the heat source temperature Tg, the 
pressure and composition of stream in the generator can be 
determined. Generally, the enthalpy of a multi-component fluid is 


h= YxXjhj+Amixh (18) 
or 
T 
h= X Cpj dT +Amixh (19) 
To 


where C,, is the heat capacity of species j of stream, Amixh is the 
mixing enthalpy of the system, which is usually ignored. To is the 
reference temperature, usually treated as 298.15 K. 


5.2. Performance assessment of absorption cycles 


Table 19 lists the published COP of refrigerant/IL systems for 
single-effect absorption cooling cycle. Deviations among the 
calculated results exist because different calculations are based 
on different property data, thermodynamic models, and operating 
conditions, e.g., heat source temperature. 


Table 19 
COP of refrigerant/IL systems for single-effect absorption cooling cycle. 


System Ref. TeC TaPC TeC THC f coP 

H20/[BMIM]BF4 [49] 10 30 40 100 13.0 0.544 
H20/[EMIM]BF4 18.2 0.525 
H20/[EMIM]EtS04 13.57 0.569 
H20/[DMIM]DMP 5.32 0.662 
H20/[BMIM]I 23.7 0.534 
H20/[BMIM]DBP 11.17 0.532 
H20/[EEIM]DEP 12.38 0.565 
H20/[EMIM]DEP 7.75 0.588 
H20/[EMIM]DMP 8.66 0.691 
H,0/[DMIM]DMP [72] 10 30 40 80 8.77 0.829 
H,0/[DMIM]Cl 9.54 0,827 
CH,O/[DMIM]DMP [66] 10 30 40 100 3.71 0.87 

TFE/[BMIM]Br [81] 10 30 40 80 439 0.807 
TFE/[BMIM]BF4 6.07 0.780 
TFE/[EMIM]BF, 5.25 0.793 
R22/[BMIM]BFe [46] 10 30 40 100 512 0.319 
R32/[BMIM]BF 7.35 0.385 
R32/[BMIM]BF, 6.41 0.330 
R134/[BMIM]BFg 438 0.348 
R134a/[BMIM]BFe 10.66 0.254 
R152a/[BMIM]BFg 13.27 0.300 
R125/[BMIM]BFg 16.49 0.128 
NHs/[BMIM]BF4 [23] 10 30 40 100 12.98 0.557 
NH3/[BMIM]BFs 17.27 0.575 
NH3/[EMIM]T£,N 24.57 0.589 
NH3/[HMIM]Cl 14.26 0.525 
NH3/[EMIM]Ac [50] 10 30 40 100 12.55 0.573 
NH3/[EMIM]EtS04 17.55 0.485 
NH3/[EMIM]SCN 12.42 0.557 


To study the feasibility of NH3/IL systems in industrial applica- 
tions, Yokozeki et al. [50] calculated the performance of absorption 
cycle using NH3/IL working pairs. In the studied imidazolium ILs, 
the cycle performance of [EMIM]Tf2N as absorbent is the best, but 
the COP is slightly lower than that of the NH3/H20 system. Given 
that the distillation separation of the NH3/[EMIMJ]Tf2N system 
does not consume too much energy, the application of NH3/ 
[EMIM]Tf2N system may be feasible. 

Results in Table 19 show that the COP of H2O/IL systems ranges 
within 0.52-0.83, which do not exceed the traditional H20/LiBr 
working pair. However, the COP of H20/[DMIM]DMP and H20/ 
[DMIMI]CI systems are close and slightly lower than that of H20/ 
LiBr. Particularly, the H20/[DMIM]DMP system overcomes the 
crystallization and corrosion problems of traditional working pair 
to some extent, and broadens the operating range. For HFC/IL 
systems, the COP ranges within 0.12-0.39, which are still relatively 
low but may be used in absorption/compression hybrid cycle. For 
NH:/IL systems, the COP ranges within 0.48-0.59, which are lower 
than that of the traditional NH3/H20 working pair, but the 
difficulty of separation is relieved. 


qc = fCp(T¢ —Tin) + Aevahtr, +Amixh (20) 


where Qc is the heat load of unit mass refrigerant produced in the 
generator, which consists of three parts. fC,(Tc-Tin) expresses the 
heat needed to warm the supercooled solution from entering the 
generator to the generation temperature, Ayal is the latent heat 
of evaporation of the refrigerant at Tc, and Amixh is the mixing 
heat of liquid refrigerant and absorbent. A larger f causes larger 
heat load of evaporator, thereby reducing the COP of the cycle. 
Fig. 8 shows that with the increase of generation temperature, 
the COP of various working pairs sharply increases initially, then 
becomes stable, and slightly declines finally. f is infinite when the 
generation temperature reaches the minimum, the required gen- 
eration heat is infinite based on Eq. (20), and the COP of the cycle is 
zero. With the increase of generation temperature f decreases, COP 
sharply increases and then smoothens. With further rise of 
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Fig. 8. COP of seven working pairs versus Tc for single-effect absorption 
cooling cycle. 
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Fig. 9. Effect of Te on COP. 


generation temperature, the temperature difference (To—Tin) 
increases, the generated heat increases, and then the COP of the 
cycle slightly decreases. At the stationary stage, the working pairs 
consisting of TFE and three different ILs, whose COP ranges 
between those of H20/LiBr and NH3/H20 working pairs, have 
wider operating range. Among these working pairs, the COP of 
the TFE/[BMIM]Br system is the highest. At the stationary stage, 
the COP of H20/[DMIM]JDMP and H20/[DMIMI]ClI are close but 
slightly lower than that of the H20/LiBr working pair. Given that 
the H,0/[DMIM]DMP system has no crystallization limits and has 
the widest operating range, the operational safety greatly 
improves. The operating range of the H2O/LiBr system comes 
second, and the operating range of the H20/[DMIMICl system is 
the narrowest. As an absorbent, [DMIM]DMP has a higher solubi- 
lity with the refrigerant, e.g., H20, and a higher boiling point than 
refrigerant. Simulation results also show that the H20/[DMIM] 
DMP system is expected to become an absorption cycle alternative 
working pair. 

For two preferred suitable H2O/IL systems H20/[DMIM]CI and 
H20/[DMIM]DMP, the effects of evaporation temperature on the 
system COP, system exergy efficiency (yg), and f are analyzed and 
compared with the traditional working pair H2O/LiBr. Fig. 9 shows 
the effect of the evaporation temperature on the COP, which is 
found to follow the order H20/LiBr > H20/[DMIM]DMP > H20/ 
[DMIM]Cl, and are all higher than 0.78. With the increase of 
evaporation temperature, the COP of the system increases, this is 
because that when the evaporation temperature is low, the 
corresponding evaporation pressure is low, f is high, and COP 
is low. 

Fig. 10 shows the effect of the evaporation temperature on ye. 
The effect of Tg on COP is opposite to that of ng. COP increases with 
the increase of Tg, whereas ng decreases with the increase of Tr. 
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Fig. 11. Effect of Tg on f. 


This is mainly because COP considers only the amount of cooling 
capacity, while the yg considers not only the amount of cooling 
capacity but also the grade of cooling capacity. In other words, the 
cooling load and cooling temperature are all important for the 
absorption cooling cycle, a lower cooling temperature results in a 
higher grade of cold, so yg decreases with the increase of Tg. Fig. 10 
shows that for the three studied working pairs, ng follows the 
order H20/LiBr > H20/[DMIM]DMP > H20/[DMIM]Cl. 

Fig. 11 shows the effect of the evaporation temperature on f, 
f values of H,0/[DMIM]CIl and H20/[DMIM]JDMP systems are 
higher than that of the H2O/LiBr system, and f decreases with 
the increase of evaporation temperature. f also markedly affects 
COP, increasing COP needs f to be reduced. To reduce f, the 
concentration difference (w,—Wp) should be increased and the 
concentration of the rich solution should be reduced. 

As summarized above, the mixture systems of refrigerant and 
IL are novel, environmental friendly absorption cycle working 
pairs. Thus, these systems are attracting increasing attentions, 
and further research and development are necessary. 


6. Conclusions 


The research frontiers of IL working pairs for absorption cycle 
in recent years were surveyed, to make the relevant work more 
rational and efficient. Based on macroscopic properties and inter- 
molecular interactions, selection methods of IL as absorbent were 
introduced. In particular, based on the author's previous studies, 
from the perspective of action mechanism between absorbent and 
refrigerant, evaluation methods of combining the UNIFAC model 
with extreme G" criterion were proposed, and the actual research 
achievements of the method were also presented. The research 
progresses of thermophysical properties of IL working pairs were 
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summarized, such as the test method, popular system, modeling of 
vapor pressure, solubility, heat capacity, and density. Based on 
thermophysical property investigations of H20, NH3, alcohol, and 
HFC systems, as well as researches of ternary systems with IL 
added to the traditional working pairs, several potential binary 
and ternary systems were put forward. Finally, working perfor- 
mances of absorption cycle with IL working pairs were calculated 
and assessed to further validate the development potential of H20, 
NHs, alcohol, and HFC working pairs. 


Nomenclature 


A; and B; parameters of Eq. (10) (dimensionless) 

An, Bn and C, parameters of Eq. (12) (dimensionless) 
a, b, and c parameters of Eq. (3) (dimensionless) 

a; and b; parameters of Eq. (11) (dimensionless) 


CFC chlorofluorocarbon (dimensionless) 
COP coefficient of performance (dimensionless) 
CG the heat capacity (kJ kg~! K7!) 


DBP dibutylphosphate (dimensionless) 


DMEU _ N,N'-dimethylethyleneurea (dimensionless) 

DMPU _ N,N’-dimethylpropyleneura (dimensionless) 

DSC differential scanning calorimetry (dimensionless) 
E181 tetraethylene glycol dimethyl ether (dimensionless) 
EOS equation of state (dimensionless) 

f circulation ratio (dimensionless) 

G excess Gibbs function (kJ kg~!) 

HC hydrocarbon (dimensionless) 

HFC hydrofluorocarbon (dimensionless) 

HFIP hexafluoroisopropanol (dimensionless) 

hi enthalpy of stream i in Fig. 7, i=1,2,...,5 (kJ kg~') 
IL ionic liquid (dimensionless) 

kı Henry's constant of the refrigerant species (kPa) 
NMP N-methyl-2-pyrrolidone (dimensionless) 

p system pressure (kPa) 

pi saturated vapor pressure of the refrigerant species (kPa) 
Q heat load (kW) 

q heat load of unit mass refrigerant (kW kg~') 

R&D research and development (dimensionless) 


R114 1,2-dichloro-1,1,2,2-tetrafluoroethane (dimensionless) 


R125 pentafluoroethane (dimensionless) 
R134 1,1,2,2-tetrafluoroethane (dimensionless) 
R134a ___1,1,1,2-tetrafluoroethane (dimensionless) 


R14 tetrafluoromethane (dimensionless) 


R143a___1,1,1-trifluoroethane (dimensionless) 

R152a ___1,1-difluoroethane (dimensionless) 

R161 fluoroethane (dimensionless) 

R23 trifluoromethane (dimensionless) 

R32 difluoromethane (dimensionless) 

R41 fluoromethane (dimensionless) 

SHX solution heat exchanger (dimensionless) 

T system temperature (K) 

To reference temperature (K) 

TFE 2,2,2-trifluoroethanol (dimensionless) 

Tin temperature of supercooled solution entering the 
generator (K) 

VLE vapor-liquid equilibrium (dimensionless) 

w mass fraction of the absorbate species (dimensionless) 

Xi liquid phase mole fraction of species i (dimensionless), 
i=1, 2, 1 presents refrigerant, 2 presents absorbant 

Vi vapor phase mole fraction of species i (dimensionless), 


i=1, 2,1 presents refrigerant, 2 presents absorbant 
[BMIM]BF, 1-butyl-3-methylimidazolium tetrafluoroborate 
(dimensionless) 


[BMIM]Br 1-butyl-3-methylimidazolium bromide 

dimensionless) 

[BMIM]Cl 1-butyl-3-methylimidazolium chloride 

dimensionless) 

[BMIM]DBP 1-butyl-3-methylimidazolium dibutylphosphate 

dimensionless) 

[BMIM]MeSO, 1-butyl-3-methylimidazolium methylsulfate 

dimensionless) 

[BMIM]PF, 1-butyl-3-methylimidazolium hexafluorophosphate 

dimensionless) 

[BMIM]Tf2N 1-butyl-3-methylimidazolium _ bis(trifluoromethyl- 

sulfonyl)imide (dimensionless) 

[BMIM]TfO 1-butyl-3-methylimidazolium trifluoromethanesulfo- 
nate (dimensionless) 

[CipMIM]Br 1-decyl-3-methylimidazolium bromide 
(dimensionless) 

[CioMIM]CI 1-decyl-3-methylimidazolium chloride 
(dimensionless) 

[C,MIM]BF, 1-alkyl-3-methylimidazolium tetrafluoroborate, 

n=1, 2,... (dimensionless) 

[C,MIMJ]Cl 1-alkyl-3-methylimidazolium chloride, n=1, 2... 

(dimensionless) 

[Cn,MIM]DMP_ 1-alkyl-3-methylimidazolium dimethylphosphate, 

n=1, 2,... (dimensionless) 

[C,MIMJEtSO, 1-alkyl-3-methylimidazolium ethylsulfate, n=1, 

2,... (dimensionless) 

[C,MIM]MeSO,_ 1-alkyl-3-methylimidazolium methylsulfate, 

n=1, 2,... (dimensionless) 

[C,MIM]PF, 1-alkyl-3-methylimidazolium hexafluorophosphate, 

n=1, 2,... (dimensionless) 

[C,MIM]Tf2N 1-alkyl-3-methylimidazolium _ bis(trifluoromethyl- 

sulfonyl)imide, n=1, 2,... (dimensionless) 

[C,MIM]TfO 1-alkyl-3-methylimidazolium trifluoromethanesul 
fonate, n=1, 2,... (dimensionless) 

[DMEA]Ac N,N-dimethylethanolammonium acetate 


dimensionless) 

[DMIM]BF, 1-methyl-3-methylimidazolium tetrafluoroborate 
dimensionless) 

[DMIM]CI 1,3-dimethylimidazolium chloride 
dimensionless) 

[DMIM]DMP 1,3-dimethylimidazolium dimethylphosphate 
dimensionless) 

[EEIM]DEP 1-ethyl-3-ethyl-imidazolium diethylphosphate 
dimensionless) 

[EMIM]BEI 1-ethyl-3-methylimidazolium bis(pentafluoroethyl 
sulfonyl)imide (dimensionless) 

[EMIM]BF, 1-ethyl-3-methylimidazolium tetrafluoroborate 
dimensionless) 

[EMIM]Br_ 1-ethyl-3-methylimidazolium bromide 
dimensionless) 

[EMIM]Cl1-ethyl-3-methylimidazolium chloride 
dimensionless) 

[EMIM]DEP 1-ethyl-3-methylimidazolium diethylphosphate 
dimensionless) 

[EMIM]DMP_ 1-ethyl-3-methylimidazolium dimethylphosphate 
dimensionless) 

[EMIMJEtSO, 1-ethyl-3-methylimidazolium ethylsulfate 
dimensionless) 

[EMIM]Tf2N 1-ethyl-3-methylimidazolium _ bis(trifluoromethyl- 
sulfonyl)imide (dimensionless) 

[EMIM]TfO 1-ethyl-3-methylimidazolium trifluoromethanesulfo- 


nate (dimensionless) 

[HMIM]BF, 1-hexyl-3-methylimidazolium tetrafluoroborate 
(dimensionless) 

[HMIMJCl 1-hexyl-3-methylimidazolium chloride 
(dimensionless) 
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[HMIM]PFs 1-hexyl-3-methylimidazolium hexafluorophosphate 


(dimensionless) 


[HMIM]Tf,N 1-hexyl-3-methylimidazolium _ bis(trifluoromethyl- 


sulfonyl )imide 
(dimensionless) 


[HMIM]TfO 1-hexyl-3-methylimidazolium _ trifluoromethanesul- 


fonate (dimensionless) 


[HydeMIM]BF, 1-(2-hydroxyethyl)-3-methylimidazolium tetra 


fluoroborate (dimensionless) 


[Ni1i20ny]Cl choline chloride (dimensionless) 
[P4444]TMPP tetrabutylphosphonium bis(2,4,4-trimethylpentyl) 


phosphinate (dimensionless) 


[P8111]TMPP  trimethyloctylphosphonium bis(2,4,4-trimethyl 


pentyl)phosphinate (dimensionless) 


Greek Letters 


Vi activity coefficient of species i (dimensionless), i=1, 2, 
1 presents refrigerant, 2 presents absorbant 

Aevalrg latent heat of evaporation of the refrigerant at Te 
(kJ kg~*) 

Amixh mixing enthalpy of the system (kJ kg~!) 

Ne system exergy efficiency (dimensionless) 

p density (g cm~?) 

yı absorption potential of the refrigerant species 
(dimensionless) 

Superscripts 

C combinatorial term (dimensionless) 

R residual term (dimensionless) 

RL Raoult's law (dimensionless) 

oo infinite dilution state (dimensionless) 

Subscripts 

A absorber (dimensionless) 

Cc condenser (dimensionless) 

E evaporator (dimensionless) 

G generator (dimensionless) 

max maximum value (dimensionless) 

p poor solution (dimensionless) 

r rich solution (dimensionless) 
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